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A canonical formulation of coupled classical-quantum dynamics is presented. The theory is 
named symmetric hybrid dynamics. It is proved that under some general conditions its predictions 
are consistent with the full quantum ones. Moreover symmetric hybrid dynamics displays a fully 
, consistent canonical structure. Namely, it is formulated over a Lie algebra of observables, time evo- 

' lution is unitary and the solution of a hybrid type Schrodinger equation. A quantization prescription 

I from classical mechanics to hybrid dynamics is presented. The quantization map is a Lie algebra 

. . isomorphism. Finally some possible applications of the theory are succinctly suggested. 
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Introduction. The interaction between the space time structure and the matter/energy content of the universe 
is described by General Relativity. The fact that all forms of matter/energy interact with the space time structure 
lies at the foundations of the theory. A pertinent question is then what is the form of the gravitational interaction 
for quantum fields In the absence of a theory of quantum gravity several attempts to describe the space time 

■ - matter/energy interaction by a semiclassical formulation of gravity have been made In such a theory the 

, space time structure is described by a classical metric while the matter content is described by standard quantum 
observables. An interesting debate has been whether these semiclassical formulations will reproduce the predictions 
of quantum gravity in the appropriate limit ||^J^ . 

This is one of several examples where what seems to be missing is a consistent general formulation of coupled 
classical-quantum dynamics. That is, a theory derived as the appropriate limit of quantum mechanics that provides 
a consistent description of a general interaction between classical and quantum subsystems |^,||J^ . 
,. It seems reasonable to expect that such a theory should satisfy the two following prerequisites: a) under some general 
' conditions - concerning the initial data and the dynamical structure of the particular system under consideration - its 
^ , predictions should be consistent with the full quantum ones, and b) the theory should display a consistent dynamical 
structure. 

^ In |p|, ^Jll| two different proposals for a theory of coupled classical-quantum dynamics were postulated and mo- 

' tivated in terms of the thus resulting properties. It was also proved ||l2| that the theory presented in [^,0 can be 
, derived as the appropriate limit of quantum mechanics. Unfortunately, the two proposals do not display 

' a consistent canonical structure. The theories are not formulated over a Lie algebra of observables and thus time 
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evolution is problematic at several levels [13ll4|. In fact, it was later proved that standard classical and quantum 
mechanics cannot be consistently coupled |15,]l6|. 

Recently a new formulation of hybrid dynamics was presented by Diosi, Gisin and Strunz The starting point 
of this proposal is quantum mechanics in the Schrodinger picture. Hybrid dynamics is then derived using coherent 
state methods under some assumptions concerning the initial data and the dynamical behaviour. However, it is not 
clear in that context whether it will display a fully developed canonical structure. 

In this letter we shall present an alternative proposal for coupled classical quantum dynamics. The starting point 
will be the full quantum formulation (in the Heisenberg picture) of a general dynamical system. Using a specific order 
dequantization jl^ we will be able to derive symmetric hybrid dynamics. The procedure automatically ensures that, 
firstly, the theory is consistent with quantum mechanics in the appropriate limit and secondly, it displays a standard 
canonical structure. 
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Canonical Structure. Let us then start by establishing the conventions and assembhng some general results. Let 
us consider a generic (iV + M)-dimensional quantum system with fundamental observables {qk,Pk) {k = 1, ■ ■ ■ , N + M) 
or sucinctly just Ok {k — 1, • • • , 2(A^ + M)). These operators act on the Hilbert space Ti. — Hi The Hilbert 

space Hi is spanned by the eigenvectors \zi, - ■ ■ ,zn > of the classical sector observables Oi {i — 1, - ■ ■ , 2N) or {qi,Pi) 
(i = 1, • • • , N). The remaining subspace H2 is spanned by the eigenvectors \wi, - ■ ■ , wm > of the quantum sector 
observables Oa (a = 2N +!,■■■, 2{N + M)) or {qa,Pa) {a ^ N + I, ■ ■ ■ , N + M). 

The algebra of linear operators of the quantum system is generated by the set: 

g = [dk, k^\,---,2{N + M)Y (1) 

This algebra is an infinite dimensional complex vector space. In p^], we saw that a possible basis for this vector space 
is given by the set of completely symmetric operators: 

B={dk,...k^^(dk,---dA ; l<ki,---,K<2{N + M)-n<^u\. (2) 



Since any observable from the classical sector commutes with all observables of the quantum sector, the set B can be 
written as the set of elements of the form: 



On-»„Oa,...„„ = (d,, •■•d,„) (dai---da„) |, (3) 



B 

with 1 < ii, • • • , in < 27V; 2N + 1 < ai, • • • am < 2{N + M); n,m ^ M. These elements are linearly independent and 
generate all elements in the algebra A{Hi ® H2). Therefore B constitutes a basis for A{Hi ® H2)- 

Let us now define the symmetric half-dequantization map. The classical phase-space associated with the classical 
sector is denoted by T*Mi. 

Definition 1: Dequantization Vg 



Vg -^- : A{Hi (g) H2) ^ 5 = {/ : T*Mi ^ ^(^2)} 



1) Vg'^' is a linear map, 

2) Vf ■«■(!) ^Icig, 

o\ yH.Q. /a a \ _ q. a 

^/ ^ S \ '''I ■■■^Ti^Ctl'-'CKm j ^ll---ln^ai---0(m ^ 

where Oi,...,^ = O^, ■ ■ ■ 0,„ e T*Mi. 

The dequantization V^ '^' satisfies the following properties: 
1) Dequantization of the product: Let A,B & A. We have: V" '^ {A ■ B) = V^ '^ {A) ® V" '^ {B), where (g) is an 
extension of the well-known ^-product [p"7|-^2| and is given by: 



SxS^S] A® B = Aexp ( ) B, (4) 
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where A = V^-^ {A)^, B = V" '^ {B) and 



N 



~ i^i I ^"p^ ] ' 

The product (g) is associative, distributive with respect to the sum and has a neutral element. Therefore the set 

(5, (g), is a complex ring. 

2) Dequantization of the hermitean conjugate: 

/ / \ t\ 

\ (6) 
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^Notice that A can be seen as an element of A{H2) of the form A = J2„ ^n^n , where yl^ G T*Mi and G A{H2)- 
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where we considered the expansion of A: A = J^^^i Iljli -^ij ^'^'^ -^ij ^ ^• 

3) Dequantization of the bracket: Using the map V^'^', we can define a new bracket structure in S. Let A,BgA. 
We have: 



7 ]]m '■ S X S ^ S; 



M 



The exphcit form of the new bracket, which is an extension of the Moyal bracket is given by: 



A,B 



M 



= A(i)B-B(g,A = Ae^l> (f J') B-Bexp (f- j") A 



A,B 



(7) 



If A, B e A(T*Mi) -the algebra of classical sector observables-, then the bracket [[, J]^^ reduces to the Moyal bracket. 
Likewise, for A,B £ A{H2), the bracket reduces to the quantum commutator. Futhermore if we disregard the terms 
of order h"^ or higher the bracket reduces to the Boucher- Traschen bracket ^ . The properties of the new bracket 
follow immediately from the properties of the full quantum one: it is antisymmetric, linear and satisfies the Jacobi 
identity. Therefore, {S, [[, ]] jv^) is a Lie algebra. Moreover, the map V^'^' is a Lie algebra isomorphism. A quantization 
prescription for symmetric hybrid dynamics is given by the symmetric half quantization map which is defined 



by A 



H.Q. 



total 



V" '^ (where V^' 



total 



: A{Hi (8) H2) — > A{T*Mi ig) T*M2) named total dequantization map, is the 



extension of the map Vg '^ to the case where the entire system is dequantized ll^). The map Ag '^' is also a Lie 
algebra isomorphism. 

4) Dynamical Structure: Consider the algebraic structure of symmetric hybrid dynamics. We can easily obtain the 
dynamical equations of the theory. They follow directly from dequantizing the full quantum equations. The time 
evolution of a general observable A is thus given by: 



A 



1 r 

ih 



A.H 



M 



(8) 



where H = 
is given by: 



Vg ''^'{H) or, alternatively, H = Ag ''^'{H) is the hybrid hamiltonian. The formal solution of this equation 



1 / it 



ri=0 



H,A 



M 



M 



(9) 
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Alternatively, the time evolution of A is generated by an unitary operator U : 

A{t) ^U(ty (g>A{0)(E)U{t), 
where U{t) is the solution of the hybrid Schrodinger equation: 



(10) 



dU ~ ~ 
in—=H(g)U, U{0) = 1cIq, 



(11) 



and satisfies W = U ^. Time evolution is a canonical transformation. More generally, all unitary transformations 
are canonical transformations: 



A,B 



®U = 



M 



>B(E)U 



M 



(12) 



and so the dynamical structure of the algebra (5, C>5) is invariant under unitary transformations. 

Predictions. Let the initial data for a general hybrid system be given by the initial wave function >S 7^2 
for the quantum sector plus a set of values {0^,5i),i = 1..2N for the classical sector (where Si are the classical 
error margins associated to the initial values 0° of the classical sector observables Oi). The aim now is to obtain 
physical predictions for the time evolution of the system. The naive procedure would be to determine a set of 
eigenvectors of Ok{t) = Ok{Oa,0^ ,t) spanning the Hilbert space 7^2 (let these eigenvectors be \bk,rnk >, so that 
Ok{t)\bk, nik >— bk\bk, mk >, where bk is the associated eigenvalue and ruk is the degeneracy index) and then assume 
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that the predictions of hybrid dynamics, for the outputs of a measurement of the observable Ok at the time t, consist 
of the set of values bk with associated probabilities p{bk) = X^m^ I < bk,mk\(j)^ > 

However there is no reason to believe that such predictions are physically valid. The reason is straightforward. The 
real physical observable is Ok (t) = Ok (Oa ,Oi,t) obtained using the full quantum formulation of the dynamical system. 
Ofe(t) has the eigenvectors \ak,nk >, where ak is the associated eigenvalue and rik is the corresponding degeneracy 
index. The physical predictions for the output of a measurement of the observable Ok are given by the values Ofc with 
probability p{ak) — J2nk I ^ Ofc,nfc|0 > P where \(f> > is the total wave function describing both the classical and the 
quantum sectors at the initial time. Clearly there is no reason why these predictions should be the same as the ones 
obtained by using the operator Ok- 

Our best chance will be to use the operator Ok to obtain some knowledge (but not the complete knowledge) about 
the outputs of a measurement of Ok ■ This can be done through a procedure similar to the one presented in [ p^ . Let 
us then summarize the main points of that method: 

i) The first step is the following: quantum and hybrid dynamics provide two alternative descriptions of the initial 
time configuration of the dynamical system. The true, physical, description is the quantum one given by a total wave 
function > that we assume, for simplicity, to be of the form {(p >~ {c/)'^ > \(f)'^ > where >g 7^2 and >e Hi- 
These two wave functions describe the quantum and the classical sector initial time configurations, respectively. On 
the other hand hybrid dynamics provides only an approximate description of the initial time configuration. Concerning 
the quantum sector, we are given exactly the same description: >, whereas for the classical sector, we only have 
the set of values (0°, Si) available. Our main task is then to understand under which conditions the two descriptions 
of the classical sector are consistent. Clearly, if they are not consistent, then we can not expect hybrid dynamics 
to yield sensible predictions. This problem was studied in detail in |Q in the context of the semiclassical limit of 
quantum mechanics and the result turned out to be a set of criteria establishing a notion of classicality. These criteria 
were then used in JT2| to test the physical validity of a proposed theory of hybrid dynamics originally suggested in 

10- 

The classicality criterion consisted of a set of conditions relating the classical and the quantum description of the 
initial time configuration of a general dynamical system that, when satisfied, ensure that the classical predictions be 
consistent (in some precise sense) with the quantum ones at all times. Here we shall also use this classicality criterion. 
We thus impose that > be first order classical (i.e. it should satisfy the first order classicality criterion) with 
respect to the classical data {0^,Si) (notice, however, that higher order classicality criteria, imposing more stringent 
conditions on the functional form of the classical sector initial data wave function, can also be used |l^,|2^). This 
means that \(f>'' > should satisfy: 

< S(^,„, 10^ >,S,J\E{S,„,\<j>'^ >,S,J >< 5%^ (13) 

where \E{Si^ , >, Si^) > is the error ket associated to the sequence of observables Si^ — Oi^....Oi^ ; 1 < ii, .., < 
2N] neJ\f, defined by: 

\E{s,^,\cj,- >,s,j >= {d,,~ol)....{d,^-olM >, (14) 

and Si^ is any sequence of observables Si^ = Oi^ ■■■■Oi^ such that: 

d^dk{t) _ d^dkit) 



^S^„ do,,.. ..do. 



^0, (15) 



for any of the observables Ok{t),k = 1..2{N + M). Moreover 6si„ — 5i,....5i^. 

In summary, a general quantum system admits a hybrid description with initial data (If/)*^ >,0^,(5i) if the initial 
data wave function \(t)^ > describing the classical sector in the full quantum description admits a proper classical 
description. This in turn means that \(f>'^ > should satisfy the set of conditions (13) for all sequences satisfying (15). 

ii) In the second step of our approach, a general relation between the operators Ok{t) = A and Ok{t) = B is 
presented and then used to establish a relation between the eigenvectors of Ok{t) and the eigenvectors of Ok{t). Let 
us then proceed along these lines. It can be proved |p^ , |l7| , that for B ~ V^''^'{A) = A: 



2N , 2JV 



BR 1 f)^ R 

^ - E 00-^0, - O,) + ^ E ^(^^ - - O,) + ... (16) 

-1 ^,j=i 



We now construct the set of eigenstates of B, {ipk, mk >= > >G Hi ® H2, which can be used to expand 

\(t> >, and derive the explicit form of the error ket of these states, in the representation of A, around the corresponding 
eigenvalue: 
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\E{A, I Vfe, ruk >, bk) >= (i - -B)| Vfe, mfc >= E?=i \E{di, \r >, Ot) >§§-\hk,mk> + 

(17) 

\E{di,djM'^ >,0.,0,) > ^05-|6fc,mfe > +... 

Using this quantity it is straightfoward to prove that A(y4, |-0fe,n^/c >,bk,p) =< E(A, {ipki'fnk >,bk)\E{A, IV'fejTOfe > 
,bk) >^''^ /(I -p)^^^ (where < p < 1 is a probabihty) satisfy: 

1 9_B^ 95 

A(i,|Vfe,mfe >,bk,p) < )i/2 Z]l <^'''^''\dada^^'"^'' ^ 1^^^^' + - 

where we exphcitly used the requirement (13). A{A,\'ipk,mk >,bk,p) is named the spread of the state \ijjk,mk > 
in the representation of A and can be used to study the properties of IV'/o^Zfe > in that representation. Namely it 
can be proved that in this representation the state \tpk,'mk > has at least a probability p confined to the interval of 
eigcnvahics of A, / = [6^. — A(A, j^/ife, m^, >, 6fe,p), 5^ + A(A, |V-'fe, "^-fc >,bk,p)]- This result is vahd for all < p < 1. The 
functional form of A suggests that we write it as A{A, \'il}k,mk >,bk,p) = <5b/(1 —pY^^ where 6b =< E{A, |^/'ft,mfe > 
,bk)\E{A, \tl)k,mk >,bk) >^l'^ is named the error of B and is of the size of a classical error margin. 

iii) Using the former relation between the eigenvectors of B and the eigenvectors of A we are able to derive a 
relation between the probabilities in the representation of B and those in the representation of A and thus to obtain 
predictions for the outputs of a measurement of A using only the knowledge about the hybrid operator B and its 
eigenvectors \bk,mk >■ Let us then consider an arbitrary interval of eigenvalues of A: Iq = [a*' — D, a^ + D], D > 2A(p) 
and let P{ak G Iq) be the probability of a measurement of the observable A yielding a value inside the interval Iq. 
The full quantum mechanical predictions are given by: 

P(afee/o)= J2 \<'P\ak,nk>\^. (19) 
Now P{ak e Iq) can be expanded using the eigenvectors of B. We have: 

2 

(20) 



P{ak e /o) = J2 



X] < <P\'^k,mk ><ipk,mk\ak,nk > 

mk,bk 



where we used the fact that \(p >= ^j,^ < Vfci'Ti/cl^ > I'^'k^'^k >■ Using the properties of the states |V'fe,'Ti/c > in 
the representation of A we can finally prove, through a quite long calculation, that: 

1 - {P{bk i IminY'^ + (1 - pY'^Y < P{ak e /o) < {P{bk e /„„.) V2 + (1 - p) V4}2 (2I) 

where Imin = [a° - {D - 2A), 0° + (D - 2A)], Imax = [a° -{D + 2 A), a° + {D + 2 A)] and A = A(i, mfe >,bk,p) 
is given by (18). Notice that we can make (1 — p) as small as desired. However, this will affect the value of A and 
thus the range of the intervals Imin and Imax ■ 

To see this explicitly let us make p = 0.99. A in (18) is given by A = lOiJs where 63 =< E{A, |V'fe,'Tifc > 
,bk)\E{A, \'iljk,m,k >,bk) >^/^ is of the order of magnitude of the classical error margins. Using this value of the 
spread we can state that, in the representation of A, the eigenvalue of B, \ipk^ "Zfe > has at least 99% of its probability 
distribution confined to the interval [6fc — A, 6fe + A]. The predictions for the outputs of a measurement of A are then: 

1 - {P{bk i Im^nf'^ + 0.32)^ < P(afc G /„) < (P(6fc £ Imaxf'^ + 0.32)^ 

(22) 

=^ P{bk e Imin) - 0.74 < P{ak e Jo) < P{bk € Imax) + 0.74 

where Imax, Imin = [do — {D ± 20Sb), ao + {D± 20^s)]- An error of 74% in a prediction of a probability is, of course, 
inadmissible. 

Let us now increase the value of p and see what happens. Let p = 0.9999. A is now given by A = IOO^b and the 
predictions for the outputs of a measurement of A: 

1 - {P{bk ^ Imin)^'^ + 0.1)2 < p(^^ g j^^ < (p(j,^ g Imax)^'^ + 0.1)^ 

(23) 

=^ P{bk e Imin) - 0.21 < P(afe e /o) < P{bk G Imax) + 0.21 
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An error of 21% is more reasonable. Notice however that the difference between the range of Iq and Imax,Imin 
increased considerably - from 206b to 2006 b- We can continue to increase the value of p but this will affect the value 
of A and thus it is not likely that we might obtain more precise predictions. 

We should have expected symmetric hybrid dynamics to provide imprecise predictions since the theory uses classical 
data to describe the initial time configuration of one of its sectors and classical data are imprecise in nature. However, 
the degree of imprecision is very large which is due to the fact that the conditions imposed on the classical sector 
initial data wave function are the least restrictive possible. It is worth pointing out that more accurate predictions are 
possible either by completely reformulating the procedure by which symmetric hybrid dynamics makes predictions for 
the outputs of a measurement of the full quantum operators Q or, albeit using the former procedure, by imposing 
more restrictive conditions on the classical sector initial data wave function . 

To see what happens in this latter case let us assume that the classical sector initial data wave function satisfies the 
2nd order classicality criterion. This means that 10"^ > should satisfy the conditions (13) for all 2nd order sequences 
of the form Si^ = (5^/^,5'^//) where Si'^ and Si'^ are any of the 1-order sequences determined in (15). In this case it 
can be proved ||l|| that A in (18) is (with a good aproximation) given by A =< E\E >^/2 /(I -p)^/^ = fc/(l -p)^^* 
and, moreover that the new spread allows for the same statement concerning the confinement of the eigenvectores of 
B. The predictions for the outputs of a measurement of A become considerably more precise: 

1 - {Pibk i hmnfl^ + (1 - vfl^l-i^'^Y < P{ak e lo) < {P{bk e Ima.y^' + (1 - pf^VS'^'V (24) 

Let us make, once again, p = 0.9999. This time A = 5^/(1 — p)^^"^ = 106b and: 

1 - iP{bk i /rmrO'/' + 0-02)2 < p(^^ ^ j^) < (p(5^ ^ /„,,)l/2 + 0.02)^ 
=^ P{bk e Irmn) ~ 0.04 < P(afe G /q) < P(5fe G /,„a.) + 0.04 



(25) 



That is an error of at most 4%, with the difference between the range of /q and /mm, imax decreasing to 206b- If we 
keep increasing the degree of the classicality conditions to be satisfied by the classical sector initial data wave function 
we will certainly increase the degree of precision of the hybrid predictions. However, we will also narrow the range of 
validity of hybrid dynamics to those systems with a classical subsystem satisfying higher order classicality criteria. 
Example. To illustrate the former results let us consider the simple hybrid system described by the hamiltonian: 

H =\^{x^ + q^) kxQ (26) 

where (g, x) are the canonical variables of a classical harmonic oscillator with m = w = 1, {Q, P) arc the fundamental 
observables of the quantum system to which the harmonic oscillator is coupled and k is the coupling constant. By 
solving the system of equations (8) we obtain the time evolution of the fundamental observables: 

q{t) ^ q{0) cost + {x{0) + kQiO)} sint 

i{t) = -g(0) sint + a;(0) cost + A:{cost - 1}Q{0) 

Q{t) = g(o) 

P{t) = P{0) - k{q{0) cost + x{0) sint} - k^{smt - t}Q{0) (27) 

where {q{Q), x{Q)) and ((5(0), P(0)) are the initial time classical and quantum sector fundamental observables, respec- 
tively. Finally the L-order spreads (18) (L e Af) are given by: 



^L{q,p) 


= <5^/(l - 


-p)V2L: 


= {\COSt\6g{0) + 


|sint|5,(0)}/(l-p)i/2i 


^l{x,p) 


- &/(l - 




= {\smt\6g{0) + 


|cost|<5,(0)}/(l-p)i/2i 


^l{Q,p) 


- V(i 




= 




Al(Ap) 


= Sp/{l 




= {\kcost\6q{0) - 


f |fcsint|4(0)}/(l-p)i/2^ 



(28) 

where < p < 1 is a probability, Sq{0), 6x{0) are the initial data classical error margins and the results are valid up 
to the order of classicality L of the classical sector initial data wave function \4>'^ >. In particular they are valid for 
the first and second order spreads L = 1,2 (if the initial data wave function > is first or second order classical) 
that were previously used to obtain the general predictions (21,24). 
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Conclusions. A canonical formulation of coupled classical and quantum dynamics was presented. The theory 
satisfies an interesting set of properties. On the one hand, and most importantly, for a general dynamical system 
satisfying some general conditions concerning its initial data and dynamical behaviour, the predictions of symmetric 
hybrid dynamics are consistent with the predictions of full quantum mechanics (and, in fact, can be used to obtain 
full quantum predictions) and are therefore physically valid. 

On the other hand, symmetric hybrid dynamics displays a fully consistent canonical structure: i) it is formulated 
over a Lie algebra of observables, ii) time evolution is unitary and probabilities are positive defined, iii) the theory 
admits a set of general canonical transformations which are generated by the action of unitary operators and iv) time 
evolution, in particular, is a canonical transformation. 

Moreover, the limit of symmetric hybrid dynamics when the classical sector does not exist is quantum mechanics. 
However, one should notice that the limit when the quantum sector does not exist is not standard classical mechanics. 
Instead this limit is symmetric classical mechanics, an alternative theory of classical mechanics that was proposed 
and studied in , and this is the reason why the results of this letter do not contradict those of . Therefore 

symmetric hybrid dynamics properly generalizes both quantum and symmetric classical mechanics. Futhermore, the 
quantization prescription from symmetric classical mechanics to symmetric hybrid dynamics is just a Lie algebra 
isomorphism and is thus not plagued with ordering ambiguities. 

Finally, notice that the limit of symmetric classical mechanics when ?i — > is just standard classical mechanics. 
Likewise, in this limit, the bracket structure and thus the dynamical structure of symmetric hybrid dynamics are 
just the ones originally proposed by Boucher and Traschen in |5|. This theory was later proved to be consistent with 
quantum mechanics under the same assumptions that were made in this letter in what concerns the initial data 
but using an approximation procedure that discarded the contributions of the terms proportional to or smaller. 
The results of the present letter corroborate, by straightforward order of magnitude considerations, the procedure 
used in JT2| and thus the validity of the Boucher- Traschen bracket as a possible dynamical structure (unfortunately 
ill behaved) for coupled classical-quantum dynamics. 
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